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Laser-based electron acceleration is attracting strong interest from the conventional accelerator 
community due to its outstanding characteristics in terms of high initial energy, low emittance and 
high beam current. Unfortunately, such beams are currently not comparable to those of conventional 
accelerators, limiting their use for the manifold applications that a traditional accelerator can have. 
Besides working on the plasma source itself, a promising approach to shape the laser-generated 
beams is coupling them with conventional accelerator elements in order to benefit from both, a 
versatile electron source and a controllable beam. 

In this paper we show that some parameters commonly used by the particle accelerator community 
must be reconsidered when dealing with laser-plasma beams. Starting from the parameters of 
laser-generated electrons which can be obtained nowadays by conventional multi hundred TW laser 
systems, we compare different conventional magnetic lattices able to capture and transport those 
GeV electrons. From this comparison we highlight some important limit of the state-of-the-art 
plasma generated electrons with respect to conventional accelerator ones. Eventually we discuss an 
application of such beams in undulators for Free Electron Lasers (FELs), which is one of the most 
demanding applications in terms of beam quality. 

PACS numbers: 



I. INTRODUCTION 

Particle accelerators have been one of the most impor- 
tant scientific instruments for the research activity of the 
last 50 years. Ever since the beginnings, an exponen- 
tial increase of the beam energies has accompanied the 
technological progress of the accelerators. However this 
tendency seems now destined to reach saturation pQ. In 
fact, the acceleration of charged particles relies on the 
use of radio-frequency longitudinal electric fields. They 
are generally produced in cavity resonators and their am- 
plitude is proportional to the energy gain of the charged 
particles. However, if the electric fields in the cavities 
reach some tens of MV/m, a breakdown of the material, 
producing a leak of electrons from the metal surface, will 
damage the resonators, preventing the particle accelera- 
tion. As a consequence, the energy gain per unit length 
in conventional accelerators is limited, and high energies 
can only be reached by additional acceleration modules, 
thus explaining the large dimensions of some state-of-the- 
art accelerators. 

A completely new approach toward particle accelera- 
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tion, which became feasible only in recent years, is based 
on the laser-plasma interaction [2]. In comparison with 
conventional accelerators, a plasma can generate and 
support electric fields up to TeV/m, thus giving the pos- 
sibility to accelerate a particle beam to energies up to 
few GeV in a distance of few cm instead of a hundred 
meters. Several experiments up to now have shown the 
capabilities of laser-plasma acceleration, producing par- 
ticle beams of high quality and mid to high energy [SI H] • 
Currently, laser-plasma accelerators can probuce beams 
of 1 GeV over 3 cm of acceleration length [4 . These laser- 
produced particle beams possess a number of outstanding 
properties, such as ultra-short pulse duration, high peak 
currents and excellent emittance at the plasma-vacuum 
interface. Given these unique beam properties, and the 
compactness of the acceleration scheme that will results 
in a strong advantage in terms of size and cost of the 
global accelerating infrastructure, the field of laser-based 
particle acceleration has attracted much attention. 

In this paper we present a design of a particle transport 
system based on conventional magnetic devices (solenoids 
and quadrupoles) that is able to capture and transport 
such novel beams. Moreover, we will show the present 
limitations of these beams, indicating the way to improve 
their quality. In fact, while currently achieved parame- 
ters do not allow for a realistic conceptual study yet, we 
find that our simulation studies can already give a useful 



guidance. 

In particular, Sec. [TT| reminds the peculiarities of laser 
generated electron by self-injection mechanism and sum- 
marizes the main characteristics of the beam used as 
the starting point for all the simulation studies. Sec- 
tion |III| presents a consistent definition of normalized 
emittance for beams with high energy spread and high 
divergence, typical not only of laser plasma sources, but 
also of wakefields accelerators. Conventional transport 
lines for such beams are critically discussed in Sec. |IVf 



VH while Sec. |VII| shows the weak points of using those 
beams in conventional undulators for the generation of 
FEL radiation. 



II. LASER-GENERATED ELECTRON BEAMS 

There are many physical processes in laser-plasma in- 
teraction that can be exploited to produce high energy 
beams. In what follows, we will focus on electron pro- 
duced in the bubble regime [5j |6] which occurs when 
the laser waist wq is smaller than the plasma wave- 
length X p and the laser parameter ao sufficiently high 
(e.g. ao ~ 2.5-3). Nevertheless our conclusions will be 
valid also for most other experimental setups, such as 
Laser Wake Field Acceleration LWFA in a mildly non 
linear regime or plasma wakefield accelerators. 

When a high power laser impinges on an under dense 
plasma, it propagates inside it and, by ponderomotive 
force, creates a volume almost completely depleted of 
electrons (i.e. the bubble). Due to the presence of 
the background positive ions (assumed to remain still), 
strong electric fields are present, whose peak longitudinal 
value scales as 



E z \V/m] « lOOVnoIcm- 3 ]^, 
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where no is the plasma density and ao the laser parame- 
ter. Such a scaling law assumes that the bubble is stable 
and matched and it requires the bubble to be roughly 
spherical. Thus the peak value for the longitudinal force 
(responsible for the electrons acceleration) has the same 
order of magnitude of its transverse component that gives 
the focusing strength for a relativistic electron [2] . There- 
fore the resulting focusing is much more intense than in 
any conventional accelerator. Some electrons are injected 
into the strongly accelerating longitudinal field at the 
bubble bottom: the details of injection depend on differ- 
ent physical processes and are thus hard to control ex- 
perimentally. In some regimes, injection will occur only 
once for a brief time [7j [8] , in some others, it will be 
subject to a non uniform trend, as the one considered 
hereafter, while in others it would happen continuously. 
There are also some experimental schemes where injec- 
tion can be controlled by a counter propagating laser [9] , 
at the price of a very low injected charge. In a recent 
paper, injection is started by tailoring the plasma longi- 
tudinal density [TU] . 



If the plasma is assumed to be cold, the transverse av- 
erage normalized electrons momentum at the injection 
into the bubble can reach a magnitude of order ao- Only 
electrons with a longitudinal component of the velocity 
greater than the bubble speed will be captured and accel- 
erated; therefore plasma electrons injected in the bubble 
are quite hot. 

Such electrons posses, as well, a transverse (thermal) 
momentum of the same order of magnitude of the lon- 
gitudinal one, that is few MeV/c in our numerical sim- 
ulations. Typical plasma accelerated beams (considered 
in this paper) reach energy of the order of the GeV and 
therefore their divergence approaches few mrads when 
they exit from plasma. On the other hand, typical elec- 
trons produced in a RF photoinjector posses a transverse 
momentum of few eV/c when extracted from the photo- 
cathode; the subsequent acceleration brings the beam to 
a few MeV energy and the divergence down to 10 -3 mrad. 

The high beam temperature and the lack of divergence 
reduction deeply influence the behavior of the beam as 
soon as it leaves the plasma channel. It is then clear 
also why increasing the energy achieved by electrons in 
the plasma would have beneficial effects on the average 
beam divergence as well. 

To summarize, the electron relativistic momentum at 
injection into the bubble is responsible of a high diver- 
gence and a large uncorrelated energy spread, as com- 
pared to bunches in conventional accelerators. Moreover, 
since the injection could occur along the whole plasma 
channel (typically few millimeters), we can also expect a 
large correlated energy spread, due to the very high field 
gradient. 

The presence of a strong focusing field implies that the 
beam is emittance dominated; its transverse distribution, 
both in charge density and in momentum, is expected to 
be well represented by a Gaussian. From the transverse 
envelope equation, after the plasma- vacuum interface, we 
can expect that the beam transverse size would increase, 
due to free diffraction, as [TT] : 



(2) 



where s is the longitudinal coordinate, e the geometric 
r.m.s. emittance, <Jo = ^(0) the initial size and <r = 
dcr x /ds(0) its first derivative. 

In order to design a transport system for laser-plasma 
electron beams, we will focus on the capture and trans- 
port of electrons generated by 100 TW class lasers, which 
are currently establishing worldwide in different labs, be- 
ing also a test facility for larger scale facilities currently 
in planning phase, such as 10 PW lasers. The global pa- 
rameters of the obtained bunch are reported in Table [T| 
which refers to the following structure: a 200 TW laser 
with wavelength A = 0.8 /im, contrast ratio of 10 10 , pulse 
duration (FWHM) r = 30 fs, and waist=15.5 /im, deliv- 
ering an intensity of I = 5 x 10 19 W/cm 2 , similar to 
what we would expect to obtain with the FLAME laser, 



currently under commissioning at the Frascati National 
Laboratories [12]. For the gas jet, we have considered the 
gas jet installed at FLAME, which produces a plasma of 
length 4.1 mm and electron density 3 x 10 18 cm -3 . These 
parameters where chosen to produce the highest beam 
energy with the available laser power [13 . 



Charge 


700 pC 


Energy 


910 MeV 


Energy spread a £ 


6.4 % 


Bunch length 


2 /mi 


Transverse size <j x 


0.5 /xm 


Transverse divergence a x > 


3 mrad 


Emittance e 


1.4xl0 -3 mm mrad 


Normalized emittance e n 


2.5 mm mrad 



TABLE I: Main beam parameters used for the design of the 
electron beam transport. The beam is assumed of circular 
spot and only the x transverse coordinate is considered. 



Regarding the beam distribution at the exit of the 
plasma (i.e. at the entrance of the transport channel), we 
have considered the one given by the 3D - PIC code AL- 
aDyn [14]. With the above mentioned laser and plasma 
characteristics, the longitudinal phase space bunch dis- 
tribution is shown in Fig. [I] Concerning the transverse 
phase space, the beam distribution appears almost Gaus- 
sian in the transverse position and divergence, as ex- 
pected for an emittance dominated beam; at the plasma 
exit, the beam is in a waist and its position and diver- 
gence are therefore uncorrelated. 
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FIG. 1: (Color online) Longitudinal phase space distribution 
obtained with the PIC code ALaDyn and used for the trans- 
port system design. 



III. HIGH DIVERGENCE AND HIGH ENERGY 
SPREAD BUNCHES 

In order to obtain a beam that can easily be trans- 
ported in a conventional FODO line, we have to match 
the initial bunch parameters at the exit of the plasma 
with those required by the FODO cell. 

Once the beam comes out from the laser-plasma in- 
teraction region, it suddenly passes from the extremely 
intense focusing fields inside the plasma to free space. A 
general strategy to control the beam should be that of 
using strong magnetic fields by means of quadrupoles or 
solenoids, as close as possible to the interaction point. 
However, from the parameters of Table [T| the Twiss beta 
function, fir = v^/a is as low as a fraction of /im, many 
order of magnitude lower than what is generally reach- 
able in conventional accelerators. Therefore, even if the 
initial transverse beam correlation is zero, it increases 
very rapidly in a drift. Figure |2ji shows the transverse 
phase space after a 1 cm drift. The figure has been ob- 
tained by using the well established macroparticle code 
TSTEP [T5], a derivative of PARMELA [16]. The trans- 
verse dimension at this point is about 60 times higher 
than the initial value and it keeps increasing at a constant 
rate; conventional magnets have not enough strength to 
counteract such a behavior and the beam control is dif- 
ficult. On top of that, at least a short drift has neces- 
sarily to be foreseen before using magnetic quadrupoles 
or solenoids for beam matching, because of the physical 
dimensions of the devices at the interaction region. 

The strong increase of the transverse beam size and 
the difficulty to control the Twiss parameters are not the 
only problems in designing a matching line. One appar- 
ently strange behavior of the bunch that turned out from 
macroparticle simulations is shown in Fig. |3j where the 
black curve displays the transverse normalized emittance 
in the 1 cm drift following the interaction region (the 
choice of 1 cm has been done just as a reference for the 
analysis of a general behavior of the bunch in the drift 
immediately following the interaction point). 

If space charge and wake field effects are neglected, 
as in this case, one should expect a constant emittance 
in a drift. This is true for the emittance e, since the 
transverse motion is uncoupled with others planes, and 
it is generally thought to be true also for the normalized 
emittance e n , since it is usually obtained, in conventional 
accelerators, as e n =< 7 > e with < 7 > the average 
relativist ic 7 factor. This way of evaluating the normal- 
ized emittance has been used for laser-plasma generated 
electron beams in international journal papers (see e.g. 
P3[8j[T7]). However, a more detailed analysis of the defi- 
nition of the normalized emittance shows that more care 
has to be taken when evaluating this quantity and that 
the bunches coming out from the laser-plasma interaction 
cannot be treated as conventional beams. 



In order to calculate the right expression for the nor- 
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FIG. 2: (Color online) Transverse phase space distribution 
after a drift of 1 cm obtained with TSTEP. 
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FIG. 3: (Color online) Increase of the transverse normalized 
emittance along a drift of 1 cm: theory and simulation results. 



malized emittance, let us start from its definition: 

e 2 n =< x 2 >< /3 2 7 V 2 > - < xPjx' > 2 (3) 

with j3 and 7 the particle relativist ic factors, and x and 
x' the transverse position and its divergence. If the cor- 
relation between the energy and transverse position is 
negligible, as in a drift without collective effects, we can 
write 



< /3 2 7 2 >< x 2 X x' 2 > 
- < M > 2 < xx' > 2 



(4) 



(5) 



(6) 



write 

2 < /?V > - < h > 2 
° E = ' 

which can be inserted in Eq. Q to give 

e 2 n = < 7 > 2 0% < x 2 >< x' 2 > + 

+ < £7 > 2 (< x 2 >< x' 2 > - < xx' > 2 ). 

Assuming relativistic electrons (/? = 1), we get 

4=<7> a W^+^)- (7) 

If the first term on the right hand side is negligible, then 
the normalized emittance is the usual < 7 > e value. For 
a conventional accelerator this is generally true: if we 
consider for example the bunch of the SPARC photoin- 
jector [18] we get that, at low energies (5 MeV), the first 
term gives a contribution to the normalized emittance a 
factor 10 3 lower than the emittance of the second term, 
and at higher energies this factor is even higher: at about 
150 MeV it is in the order of 10 5 . 

If we consider the laser-plasma beam of Table II we find 
that the first term on the right hand side of Eq. (71, at the 
plasma- vacuum interface, has the same order of magni- 
tude as that of conventional accelerators at low energies, 
i.e. about a factor 10 3 lower than the geometric emit- 
tance; anyway, after only 1 cm of drift, due to the rapid 
increase of the bunch size, the first term becomes pre- 
dominant, being more than three times the second one. 
Therefore the normalized emittance increases in a drift. 

At the plasma- vacuum interface, the beam is in a waist, 
thus one can write in Eq. ^ o-' =0 and e 2 /a 2 , = a 2 ,. In 
the limit s ^> 1, Eq. Q reads <j x (s) « <j x 'S and Eq. ^ 
becomes 



^<7> (cr E cr x ,s 



(8) 



The definition of relative energy spread ge allows us to 



This last equation is represented in Fig. [3] by the red 
curve and it shows a very good agreement with the re- 
sults of the multiparticle simulation code. Notice that 
for protons and heavier ionic species, the assumption of 
relativistic motion done in deriving Eq. ^ is usually not 
valid; this worsen the effect on normalized emittance and 
produces a lengthening of the produced bunch. 

The physical reason behind such a dramatic emittance 
dilution is readily understood when we realize that the 
betatron frequency of a beam critically depends on its 
energy. During the drift, each chromatic component ro- 
tates with its own velocity in the transverse phase space, 
spreading out the area occupied by the whole bunch, as is 
clearly shown by Fig. [2Jd; the resulting 'projected normal- 
ized emittance then becomes a function of both position 
and energy spread. In Fig. [4^i we can see a spectrum 
of the bunch where a few (arbitrarily chosen) energetic 
components have been outlined using different colors: the 
chromatic effect is then pictorially shown in Fig [4}d where 
the transverse phase space is reproduced after removing 
the x — p x correlation. Notice that, if the divergence x r is 
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FIG. 5: (Color online) Normalized emittance along a 
quadrupole magnetic transport line: comparison between sim- 
ulations and theory. 
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FIG. 4: (Color online) Bunch spectrum (up) and de-correlated 
transverse phase space (down). Here the free drift is about 1 
m long, in order to enhance the effects of chromatism. 



defined as x' « p x / < p > (see, for example, [19]), such 
an emittance dilution is not detected, as seen in Fig. [2^l. 

Equation ^ can be also used to predict the normal- 
ized emittance behavior in presence of magnetic fields, 
for examples quadrupoles, when there is some coupling 
between energy and transverse coordinate. In Fig. [5] 
we show a comparison between Eq. ([7]), the expression 
< 7 > e and the normalized emittance obtained with the 
macroparticle simulation code. The agreement between 
theory and simulation is worse than in the case of a sim- 
ple drift; nevertheless Eq. Q remains a valid analytical 
tool to quickly determine the behavior of the normalized 
emittance in a transport line without running macropar- 
ticle tracking codes. 

We must underline that chromatic effects could be dif- 
ferent for slice normalized emittance, which is the criti- 
cal parameter for some applications; for example in FELs 
the longitudinal slices are of the order of the cooperation 
length. If such slices are sufficiently thin, the contribu- 
tion of the correlated energy spread to the slice emittance 
is usually negligible with respect to uncorrelated energy 
spread (which depends on thermal emittance). There- 
fore, concerning the slice parameters, laser-plasma beams 
may not differ much from beams from conventional ac- 
celerators. In the upper plot of Fig. [6| a comparison 
between slice emittance at s = and s = 1 cm is shown, 



while in the lower plot we show the slice energy spread. 
It is readily apparent how slice emittance is almost con- 
served in the longitudinal range between 2 and 4 /im, i.e. 
where the slice energy spread is relatively low (< 1%). 
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FIG. 6: (Color online) Comparison of slice emittance before 
and after free propagation (upper plot) together with slice 
energy spread (lower plot). The green curve reports the slice 
current in order to show the longitudinal charge distribution, 
The number of slices is arbitrarily set to 200. 

The results here presented demonstrate that the laser- 
plasma electron beams have a behavior quite different 
from beams of conventional accelerators. In a drift im- 
mediately following the laser-plasma interaction point, 
we find the uncommon situation in which the normalized 



emittance not only increases, but with a rate as high as of 
1000 mm mrad/m. This is a rather new result that must 
be carefully considered in the design of a beam transport 
system. It depends on the combination of the RMS val- 
ues of energy spread, beam divergence and beam size. To 
match such a beam to a transport line, without spoiling 
the bunch normalized emittance, requires the optical el- 
ements to have a focal length of the order of the bunch /3 
function, which is impossible to achieve with conventional 
technology. Moreover, since the bunch high chromatic- 
ity, any possible viable matching would be valid only for 
that portion of charge possessing the proper energy: in 
fact, an achromatic focusing device is not realizable by 
magnetic quadrupoles or solenoids [20 . 

In order to overcome these problems, it is necessary 
to act on the plasma-laser source itself, because conven- 
tional magnetic transport line are inadequate, as we will 
see in the next sections. Plasma lenses could be employed 
as well [2Tj . 

IV. QUADRUPOLE MATCHING LINE 

We have investigated two ways for controlling and 
matching the beam to a FODO transport line: by using 
quadrupoles or solenoids. A quadrupole matching line is 
more commonly used in conventional particle accelera- 
tors due to the lower cost as compared to a solenoid line. 
For the design of our matching system, we have first used 
the code TRACE 3-D [22J, an interactive beam dynam- 
ics program that calculates the envelopes of a bunched 
beam, including linear space charge forces, through a user 
defined transport system. Once the beam line has been 
optimized, the magnet parameters have been inserted in 
the tracking code TSTEP to account for the beam dis- 
tribution function. 

In Fig. [7] we show the optimized beam line obtained 
with TRACE 3D. Starting from very low Twiss fix val- 
ues, of the order of 0.2 mm, by using a matching system 
of three quadrupoles, with a total length of 174 mm, 
we manage to obtain at the end of the transport line a 
bunch with Twiss parameters ar=0 and /3t=45 m. Such 
final beam parameters can be transported with a peri- 
odic FODO system without particular problems. Also 
the transverse dimensions of the bunch allow a conven- 
tional transport after the matching line. 
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FIG. 7: (Color online) Matching transport line with 
quadrupoles given by TRACE 3-D. Axis limits for the dif- 
ferent shown parameters are indicated by the normalized or- 
dinate in figure. 



For a deeper investigation of the beam dynamics, we 
have introduced the magnets parameters and dimensions 
into the code TSTEP. The bunch remains confined be- 
low 0.4 mm and the behavior agrees with the results of 
TRACE 3D. 

Even if the quadrupole matching line is quite attrac- 
tive, because of its simplicity and cheapness, two prob- 
lems arise by using this kind of transport system: one is 
a pure optical problem, the other is related to the defini- 
tion of normalized emittance. 

From the point of view of the optics, the strength of the 
quadrupole gradient necessary to contain the divergence 
of such beams reaches values as high as 5000 T/m, a fac- 
tor 10 higher than the maximum gradient currently avail- 
able and obtained using permanent magnet quadrupoles 
[23] . We tried to reduce the quadrupole gradient to val- 
ues of 500 T/m, that is the actual state of the art, and 
increase the number of quadrupoles, but we did not man- 
age to control the beam explosion. 

The second problem is related to the normalized emit- 
tance e n . Using the quadrupole matching line that has 
been optimized for the optics, even if we allow the pres- 
ence of permanent magnet quadrupoles with gradients of 
5000 T/m, we obtain that the normalized emittance at 
the exit of the line reaches values up to 250 mm mrad, 
excessively high for an accelerator. In Fig. [8] we show the 
results of the simulations obtained with TSTEP starting 
from the distribution function of Fig. [I] (blue line) and 
with an ideal Gaussian distribution (green line). As ex- 
pected, the main problem is not related to the kind of 
distribution, but to the initial beam parameters. 




FIG. 8: (Color online) Transverse normalized emittance along 
the quadrupole line obtained with TSTEP. 



V. SOLENOID MATCHING LINE 

An alternative to quadrupoles is the use of solenoids 
that, acting over a longer distance than quadrupoles, al- 
low a smoother control of the bunch transverse size. As 
in the previous section, we have first optimized the beam 
line by using TRACE 3-D. A solenoid of about 20 cm, 
located very close to the laser-plasma interaction point, 
is able to control the beam size and the Twiss beta func- 
tion. At the exit of the solenoid we obtain a beta function 
of about 88 m in both directions, and a transversely un- 
corrected bunch (olt = 0). 



We have then performed simulations with TSTEP, 
obtaining a transverse size behavior similar to that of 
quadrupoles, remaining confined below 0.4 mm. 

A transport line with solenoids is more expensive and 
complicated with respect to one with quadrupoles. How- 
ever, in this case, the intensity of magnetic field inside 
the solenoid, even though extremely high, seems feasi- 
ble with the actual state of the art technology. In fact, 
the optimization with TRACE 3-D requires a longitudi- 
nal magnetic field of about 45-50 T, an unusual value for 
conventional accelerators, since e.g. for SPARC [24] and 
LCLS [25] photoinjectors it is below 0.5 T. Such high 
field solenoids have been under investigation by using 
high temperature superconductors, and a conceptual de- 
sign of a solenoid of 45 T has been recently presented 
[26] . However, as expected, the normalized beam emit- 
tance endures the same effect of uncontrollable increase 
as the one shown in Fig. [8] At the end of the beam- 
line it reaches a value of about 200 mm mrad, not much 
different from the result with quadrupoles. 

Even though a transport system by using a strong 
solenoidal field seems feasible from the optics point of 
view, the problem of an uncontrollable increase of the 
normalized emittance still affects these laser-plasma gen- 
erated bunches due to the combination of high energy 
spread and divergence, and conventional magnet trans- 
port systems are not able to counteract this effect. We 
can conclude that to contrast the phenomenon it is nec- 
essary to act directly on the laser-plasma source, by im- 
proving the beam qualities, for example, with a laser- 
beam shaping. 



and transverse phase space. In fact, in Fig. [9] we show 
the transverse phase space after the cut with the blue 
dots representing all particles whose energy differs less 
than 2% from the average and with the red dots all other 
particles. 




FIG. 9: (Color online) Transverse beam distribution obtained 
with TSTEP after the cut with the rectangular aperture: blue 
dots represent particles having energy spread less than 2%, 
and red dots all the others. 

Since in a drift all the particles having a divergence 
greater than zero move in the right direction of the figure, 
while particles with negative divergence move toward left, 
the bunch rotates, enlarging the phase space: a cut at the 
correct position then allows to reduce the energy spread. 
The same phase space after a drift of 80 cm (i.e. at 140 
cm from the laser-plasma interaction point) is shown in 



Fig. [10] If we place here a rectangular aperture with the 
proper dimensions, we are able to cut particles having 
high energy spread. 



VI. CONVENTIONAL STRATEGIES TO 
IMPROVE BEAM TRANSPORT 

If the total beam charge is of the order of magnitude as 
in Table [T] the current is about 100 kA, a factor between 
100 and 1000 times higher than that required in pho- 
toinjectors for driving a FEL. The idea is then cutting 
the bunch, thus reducing the beam divergence and the 
energy spread at the cost of loosing some charge. The 
simplest way to do that is by using slits or apertures, 
placed in optimized points of the matching line, which 
cut transversely the beam. 

A first rectangular aperture has been positioned at 60 
cm from the laser-plasma interaction point, where the 
longitudinal magnetic field is zero. The aperture cuts 
the bunch horizontally at ±0.4 mm and vertically at ±0.3 
mm. In this way we pass from the initial charge of 700 pC 
to about 230 pC. The initial current is reduced to about 
40 kA, a still very high value, and the positive aspect is 
that the normalized emittance is drastically pulled down 
of a factor higher than 7. 

With this cut alone, we cannot affect much the en- 
ergy spread that remains almost unchanged, passing from 
6.4% to 6%. It is interesting to note however that in this 
position there is a profitable correlation between energy 




FIG. 10: (Color online) Transverse beam distribution ob- 
tained with TSTEP at 140 cm from the laser-plasma interac- 
tion point, with blue dots representing particles having energy 
spread less than 2%, and red dots all the others. 

This aperture has been chosen to be 80 fim large and 60 
fim high, allowing us to obtain a bunch with an energy 
spread of 2.74%, a charge of about 70 pC, and a peak 
current of 13.5 kA. Such a bunch has characteristics close 
to those of conventional accelerators and is not so hard 
to match to a FODO periodic line, since the Twiss beta- 
functions are about 20 cm. This matching line gives a 
normalized emittance of about 10 mm mrad, a factor 
40 lower than that obtained without cuts. Moreover, 
the transverse cuts are quite easy to realize with simple 
rectangular apertures and, despite the particle loss, we 
end up with a quite high current. 



However the energy spread can not be completely con- 
trolled. With transverse cuts, this quantity cannot be 
reduced to values less than 1% - 2% without a big loss of 
current. If we want to keep a quite high current, other 
mechanisms should be devised, e.g. a transport through 
a bending magnet or an isochronous system. In fact, Fig. 
|TT[ shows the longitudinal phase space particle distribu- 
tion after the second aperture; particles with high energy 
(> 1 GeV), far from the core of the bunch, give a high 
contribution to the energy spread, but a small one to the 
whole beam current. 
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FIG. 11: (Color online) Longitudinal phase space particle dis- 
tribution obtained with TSTEP after a second beam cut with 
transverse rectangular aperture. 

Particles with such an energy distribution could not 
be cut by using simply transverse apertures, without 
first making the beam pass through a dispersion region, 
realized with bending magnets that produce energy de- 
pendent transverse trajectories. These can eventually be 
exploited to create a transverse beam-energy correlation 
that allows to properly cut the energy distribution. How- 
ever the high beam qualities are ruined as the beam en- 
ters in a bending magnet, making again hard the control 
of the bunch characteristics. Moreover, given the high 
current, coherent synchrotron radiation can generate mi- 
crobunching. 

The procedure described above could fail to achieve 
the desired beam parameters on most shots because of 
intrinsic limits in the stability and repeatability of self- 
injected beams. In particular, laser pointing instabil- 
ity and inhomogeneities in the plasma density produce 
a spread in the beam centroid direction, with respect to 
the design trajectory, that is in the order of some mrads. 
Since the second aperture on our transport line is smaller 
than the position spread expected at 140 cm, most shots 
would just be completely stopped. An alternative strat- 
egy could be tuning the transport line on a well definite 
energy, e.g. the neatly defined energy peak seen on the 
spectrum in Fig. [4] We expect that the electrons with 
higher or lower energies, being unmatched, would be lost 
to the pipe walls, managing to produce a beam with low 



energy spread and not significantly affected by chromatic 
effects. However, this scheme can also fail, due to the 
same limits on stability and repeatability, since the neat 
peak energy inside the bunch, its central energy and the 
average energy of the whole bunch are all strongly af- 
fected by jitters typical of laser-plasma interaction. 

Therefore, we can conclude that a great amount of 
work needs to be done on the production side of the 
plasma generated beams before reaching a quality com- 
patible with conventional accelerator transport line. Ac- 
cording to the discussion above, the plasma beam qual- 
ity must be increased looking for strategies to reduce the 
energy spread, the beam divergence and beam centroid 
spread inside the plasma itself. 



VII. APPLICATION TO FREE ELECTRON 
LASERS 

In order to asses the potentiality of a laser-plasma gen- 
erated beam, we have numerically studied the propaga- 
tion of such a beam into undulators for the generation of 
hard X-rays; in the following, we deal with the beam after 



the matching discussed in Sec.|VI| The simulations of the 
output characteristics of the SASE FEL radiation, such 
as temporal and spectral profiles, have been performed 
by means of the time dependent three-dimensional simu- 
lation code GENESIS 1.3 [27 , accounting for slippage, 
diffraction, emittance and energy spread effects. The 
SPARX Hard X-ray undulator beamline (undulator pe- 
riod and RMS parameter, A n =1.5 cm and i^ J RMS'=0.907, 
respectively) [28 has been used to produce SASE FEL 
radiation at 4.6 nm at the fundamental wavelength, with 
electron beam average parameters reported in Table |ll| 
and beam current profile of Fig. [12] . 



Electron charge Q 


72 pC 


Beam energy 


881 MeV 


Energy spread (rms) 


2.74 % 


Bunch length (rms) a T 


1.6 /im 


Average current Q / {s/2ti(t t ) 


5.4 kA 


cr x , cr y 


33, 27 fim 


® x' i 


0.2 mrad 




12.8 mm mrad 



TABLE II: Average beam parameters used for the SASE FEL 
simulation. 

The saturation power along the undulator is shown in 
Fig. [13] for a bunch with Gaussian current distribution 
and with the parameters reported in Table [il] Figure [l4j 
instead, is the saturation power when the actual bunch 
distribution obtained by a TSTEP simulation is injected 
into the undulator. For the Gaussian bunch, the satu- 
ration length is much longer than the undulator length 
while, for actual one, the onset of saturation is evident 
at 15 m (being the 3D gain length about 0.7 m) because 





FIG. 12: (Color online) Beam current profile at the entrance 
of the undulator sections. 



FIG. 14: (Color online) Maximum power in the radiation 
pulse along the undulator for the TSTEP-generated beam 
distribution. 



of the slices with higher current. In the simulation we 
noticed that the peak current decreases along the undu- 
lator because of particle losses due to both a non well 
optimized matching at the undulator entrance and large 
beam divergence. Therefore the current of the bunch ac- 
tually contributing to the SASE FEL process is much 
smaller than the one shown in Fig. [12] in fact, we would 
have expected an even smaller saturation length if the 
current distribution remained constant along the undu- 
lator and equal to Fig. [12] The issue of such beams is 
that the excellent properties (e.g. brilliance) can not be 
exploited because of the difficulties in matching (due to 
the big energy spread and divergence). 



-After 15 m 
-After 24 m 

-After 31 .050 m (undulator exit) 
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FIG. 13: (Color online) Maximum power in the radiation 
pulse along the undulator for the gaussian bunch case. 

The spectral analysis for the 4.6 nm radiation is shown 
in Fig. [l5] at the saturation (15 m, black curve), after 24 
m (red curve) and at the end of the undulator (green 
curve). The emission is far from being monochromatic 
and it is centered around the fundamental wavelength, 
with a 3% radiation bandwidth spread. 

Our simulations have demonstrated that in spite of 
large emittance and huge energy spread, thanks to the 
high peak current, SASE FEL saturation occurs even in 
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Wavelength (nm) 

FIG. 15: (Color online) Spectral structure for 4.6 nm radi- 
ation from the SPARX Hard X-ray undulator beamline at 
saturation (31 m) in case of the external distribution from 
TSTEP. 



moderate undulator lengths (despite the significant por- 
tion of charge lost in the propagation in the undulator). 
However, the radiation spectrum has not the shape typi- 
cal of conventional SASE FELs (spike radiation) because 
the energy spread is too high (even higher than the FEL 
p parameter). 



VIII. CONCLUSIONS 

In this paper, we have shown that laser-plasma gener- 
ated beams possess peculiar characteristics with respect 
to the conventional accelerator beams: for example in a 
drift space, the normalized projected emittance increases 
due to the high energy spread and divergence. We have 
demonstrated that capture of such beams with conven- 
tional accelerator devices is in principle feasible: one can 
use solenoids and rectangular apertures to properly shape 
the beam, even despite some beam charge losses, and de- 
vise an isochronous beam line to further reduce the en- 
ergy spread. Applications, such as FELs (i.e. requiring 



good slice beam properties), can profit of the fact that 
slice emittance may not be diluted too much since slice 
energy spread can be quite low for such beams. Anyway 
we believe that the qualities of such beams (in terms of 
energy spread and beam divergence) must be improved 
first inside the plasma itself (e.g. properly shaping the 
plasma channel, control over injection, etc..) before al- 
lowing the use in conventional accelerator lines. 
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